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Q\ • We present an evolutionary scenario for the spin-up and evolution of binary 

millisecond pulsars, according to which the companion of the pulsar PSR J1740- 

CN ■ 5340 , recently discovered as a binary with orbital period of 32.5 hrs in the 

^ ■ Globular Cluster NGC 6397, is presently in a phase of 'radio-ejection' mass loss 

from the system. The optical counterpart is a star as luminous as the cluster 

£T) ■ turnoff stars, but with a lower T e g (a larger radius) which we model with a 

star of initial mass compatible with the masses evolving in the cluster (~ 0.85 
M Q ). This star has suffered Roche lobe overflow while evolving off the main 
■ sequence, spinning up the neutron star to the present period of 3.65 ms. At 

present, Roche lobe overflow due to the nuclear evolution of the pulsar companion 
2 ■ and to systemic angular momentum losses by magnetic braking is still going on, 

but accretion is inhibited by the momentum exerted by the radiation of the 
pulsar on the matter at the inner Lagrangian point. The presence of this matter 



around the system is consistent with the long lasting irregular radio eclipses seen 
in the system. Roche lobe deformation of the mass losing component is also 
necessary to be compatible with the optical light curve. The "radio-ejection" 
phase, which had been recently postulated by us to deal with the problem of the 
lack of submillisecond pulsars, can be initiated only if the system is subject to 
intermittency in the mass transfer during the spin-up phase. In fact, only when 
the system is detached the pulsar radio emission is not quenched, and may be 
able to prevent further mass accretion due to the action of the pulsar pressure at 
the inner Lagrangian point. We propose and discuss that a plausible reason why 
the system is expected to detach is the irradiation of the mass losing component 
by the X-ray emission powered during the accretion phase. We finally discuss 
the consequences of the binary evolution leading to PSR J1740-5340 , and its 
relation with other possible optical counterparts. 
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1. Introduction 

The widely accepted scenario for the formation of a millisecond pulsar is the recycling 
of an old neutron star (hereafter NS) by a spin-up process driven by accretion of matter and 
angular momentum from a Keplerian disc, fueled via Roche lobe overflow of a binary late 
type companion (see Bhattacharya & van den Heuvel 1991 for a review). If the NS has a 
magnetic dipole moment (typical values are /i ~ 10 26 — 10 27 G cm 3 ) the disc is truncated at 
the magnetosphere, where the disc pressure is balanced by the magnetic pressure exerted by 
the NS magnetic field (Pmag oc /i 2 r~ 6 ), where r is the generic radial distance from the NS 
center. Once the accretion and spin-up process ends, the NS is visible as a millisecond radio 
pulsar (hereafter MSP). Indeed, a common requirement of all the models of the emission 
mechanism from a rotating magnetic dipole is that the space surrounding the NS is free of 
matter up to the light cylinder radius Plc (at which the speed of a body rigidly rotating 
with the NS equals the speed of light) . 

An interesting evolutionary phase can occur during the accretion onto the NS if the mass 
transfer rate drops below the level required to allow the expansion of the magnetosphere 
beyond Plc switching-on the emission from the rotating magnetic dipole (e.g. Illarionov 
and Sunyaev, 1975; Ruderman, Shaham & Tavani 1989; Stella et al., 1994). In a standard 
Shakura Sunyaev accretion disk, it is common to identify four zones based on the different 
roles of opacities (free-free or electron scattering) and pressures (radiation or gas pressure). 
The radius of separation between the zone (B), in which the electron scattering opacity 
and gas pressure dominate, and the most external zone (C), where the free- free opacity 
and gas pressure dominate, is P BC = 4.18 x 10 7 L% 3 R 2 6 /3 m x / 3 cm (see, e.g. Burderi, King 
& Szuszkiewicz 1998), where L 37 is the luminosity in units of 10 37 erg s _1 , Rq is the NS 
radius in units of 10 6 cm, and m is the NS mass in solar masses. The typical distance 
of the inner Lagrangian point from the NS is larger than ~ 10 11 cm for orbital periods 
longer than ~ 30 hr. Close to the Roche-lobe the disc is then in zone C, where opacity is 
dominated by free-free processes and pressure is dominated by the gas contribution. In this 
case, Pdisc oc M 17//20 r~ 21 / 8 , where M is the accretion rate. The radial dependence of the 
pulsar radiation pressure is Prad oc /i 2 P~ in r~ 2 (where P sp in is the NS spin period), and 
then it is flatter than the radial dependence of the disc pressure close to the inner Lagrangian 
point. Therefore, Prad dominates over Pdisc for any t > Plc and the whole accretion disc 
is swept away up to the inner Lagrangian point L\ by the radiation pressure of the pulsar. 
During this "radio ejection" phase, the mechanism that drives mass overflow from Li can 
well be active, but the pulsar radiation pressure at L\ prevents mass accretion onto the NS. 

The occurrence of a radio-ejection phase in the pre-MSP evolution has been recently 
proposed by Burderi et al. (2001a,b), to deal with the problem of the apparent lack of sub- 
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millisecond pulsars, namely radio pulsars with spin periods below 1 ms, (Burderi & D'Amico, 

1997) and with the problem that MSP do never seem to harbor very massive NSs (see e.g. 
Taam, King & Ritter, 2000, Tauris & Savonije 1999), as it would be the case if conservative 
or quasi-conservative evolution has taken place. 

A typical progenitor system of binary MSPs is made by a NS (hereafter the primary) 
accreting from a giant, or from a star on its way to become a giant (hereafter the secondary). 
This is testified by the many MSPs with low mass (M ~ 0.2 — 0.4 M ) white dwarf (WD) 
binary companions (e.g. Rappaport et al. 1995, van Kerkwijk et al. 1996, Burderi et al. 

1998) , the remnants predicted by this type of semidetached evolution. The NS spin-up 
occurs via the accretion of angular momentum of the mass overflowing the Roche lobe. This 
mass transfer is dictated by the nuclear evolution of the secondary, and the losses of systemic 
angular momentum also contribute to determine the mass transfer rate. It is well possible 
that mass transfer during the spin-up phase suffers instabilities (see below) and the system 
temporarily detaches, allowing the pulsar to switch-on. However the companion evolution, 
which leads to radius expansion, will lead again to overflow. Mass exchange to the NS and 
spin-up will go on, until the pulsar has been so much spun up that its radiation pressure at 
the inner Lagrangian point is high enough to prevent mass accretion. Modeling of the binary 
evolution at the time of action of this radio ejection mechanism is non trivial, and an adequate 
study of the physics involved and numerical simulations are needed to adequately investigate 
this problem, as it has been done for the case of "evaporation" of low mass companions of 
MSPs (Tavani & Brookshaw 1991, Banit & Shaham 1992). The destiny of the system is 
crucially dependent on global mass radius exponent of the mass losing component, on the 
mass ratio and on the specific angular momentum of the ejected matter / (see Ritter, 1995, 
for general discussion of the stability of mass transfer). We can foresee cases in which / is 
very large with respect to the specific angular momentum of the system, and a dynamical 
phase of mass transfer can occur (Burderi et al., 2002, in preparation). Here, we will follow 
only the case in which / is the specific angular momentum at the inner lagrangian point, 
and the evolution is stationary: the orbit widens and the companion placidly expands on its 
nuclear timescale, keeping its contact with the Roche lobe. 

The main consequence of the occurrence of the radio-ejection phase, will be that the 
matter lost from the secondary will not be transferred to the NS. This phenomenon may look 
like the so called "ablation" of the companion irradiated by the wind (composed of electro- 
magnetic radiation and relativistic particles) of a MSP (Tavani 1991). The key difference 
is that, in the ablation mechanism, the process takes place once the companion is detached 
(and is evolved into a low-mass white dwarf), while radio-ejection takes place during the 
active phase of Roche lobe overflow by the secondary, and the energy requirements are much 
easier to be met. Moreover, the energy available to ablate the companion is different, being 
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the orbital binding energy in the case of radio-ejection and the pulsar rotational energy in 
the case of ablation. Observationally, in both cases we expect to be in the presence of a 
radio MSP which from time to time is obscured by the matter floating around the system. 

In this paper we propose that this postulated radio ejection phase has now been detected 
in the system containing the eclipsing millisecond pulsar PSR J1740-5340 , discovered by 
D'Amico et al. (2001a) in the globular cluster NGC6397, and identified (Ferraro et al. 2001) 
with a slightly evolved turnoff star in the sample studied by Taylor et al. (2001). 

We will provide convincing evidence that PSR J1740-5340 is an example of a system 
in the radio-ejection phase, by modeling the evolution of the possible binary system pro- 
genitor. We also discuss why we should expect that this type of system may indeed suffer 
the detachments required in order to initiate the radio-ejection phase: in fact, the X-ray 
illumination of the secondary star during the accretion phases, which leads to the NS spin 
up, is such that the mass transfer rate is probably very unsteady. When: 

i) as a result of the accretion of matter and angular momentum, the NS spin period is so 
short that, potentially, the radiation pressure of a pulsar phase would be high enough to 
overcome the pressure of the matter overflowing the Roche lobe, and 

ii) the oscillations in M are large enough to allow the millisecond pulsar to switch-on, 
the radio-ejection phase begins, leading to the appearance of the system as it looks now. 

We also show that the most luminous 'hot' objects in NGC 6397 identified by Taylor 
et al. (2001) as low mass helium WDs, indeed follow the evolution of binaries which may 
account for the evolutionary status of PSR J1740-5340 , while the less luminous stars in this 
sample might have radii too large to be helium WDs. 

2. PSR J1740-5340 

The millisecond radio pulsar PSR J1740-5340 was discovered in the globular cluster NGC 
6397 (D'Amico et al. 2001a). It has the longest orbital period (P or b — 32.5 hrs) and the most 
massive minimum companion mass (0.18 M ) among the 10 eclipsing pulsars detected up to 
now. The spin period (P S pm — 3.65 x 10~ 3 s) and its derivative (P S pm = 1.59 x 10~ 19 ), recently 
derived by D'Amico et al. (2001b), allow the determination of the NS magnetic moment 
fi 2 Q — 7.7. Because of its position with respect to the cluster center, the contamination 
of -Pspin because of the NS acceleration in the gravitational field of the cluster is negligible 
(D'Amico et al. 2001b), implying that the estimate of the NS magnetic moment is reliable. 

Ferraro et al. (2001) have extracted, from archive HST data, a star whose position 
is consistent with the position of PSR J1740-5340 (derived from one year of pulse timing 
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by D'Amico et al. 2001b). This optical counterpart shows light modulation at the same 
orbital period as the radio data. D'Amico et al. (2001b) have studied the radio eclipses 
of PSR J1740-5340 , which last for about 40% of the orbital phase at 1.4 GHz. Out of 
eclipse the pulsar signal at 1.4 GHz shows significant excess propagation delays (up to ~ 8 
ms) and strong intensity variations. Similar variations have been observed, at 400 MHz, in 
the eclipsing pulsar PSR B1957+20 close to the eclipse ingress and egress (Fruchter et al. 
1990) and in PSR J2051-0827, at 1.4 GHz, when the pulsed signal is occasionally detected 
during the eclipse (Stappers et al. 1996). These similarities suggest that in PSR J1740- 
5340 the signal is propagating through a dense material surrounding the system. In order to 
investigate this possibility, D'Amico et al. (2001b) have fitted the excess delays measured in 
two adjacent bands of 128 MHz each at 1.4 GHz. They found that the excess delays At can 
be well fitted with the equation At oc //- 2 - 02±0 - 30 that strongly supports the hypothesis that 
the responsible mechanism is dispersion in a ionized medium. In this case the corresponding 
electron column density variations are An c ~ 8 x 10 17 At_3 cm" 2 , where At_3 is the delay 
at 1.4 GHz in ms. For At_3 ~ 8 the estimated electron column density is ~ 6.4 x 10 18 cm" 2 . 

The eclipsing radius is Re ~ asin(0.47r), as the eclipse lasts for ~40% of the orbit 
(D'Amico et al., 2001b). It is 4.4 x 10 11 cm, if we take m 1 = 1.8M for the NS mass 
and m 2 = 0.45M Q for the secondary mass (see below) and P orb ~ 32.5 hr (a is the orbital 
separation). This radius is larger than the Roche lobe radius of the secondary (~ 1.3 x 10 11 
cm). This means that the eclipsing matter is beyond the gravitational influence of the 
companion star and must be continuously replenished. We can estimate a rough order of 
magnitude of the necessary mass loss rate from the secondary, M, by assuming spherical 
symmetry (which is not consistent with the randomly variable signal intensity shown by the 
radio data). We further take the minimum electron column density which, for a system seen 
edge-on, occurs when the NS is in front of the secondary, so that: 



where 7 is the fraction of ionized gas, -Rrli is the radius of the primary Roche lobe, 
X ~ 0.7 and Y ~ 0.3 are the hydrogen and helium mass fractions respectively, and m p is 
the proton mass. Using for the density p the continuity equation: 




(1) 



A7rr 2 pv = M 

where f W md is the speed of the outflowing matter, a simple integration gives 



(2) 



M_ 10 < 1.17 x n 20 1~\X + 0.5F)- 1 v 8 (mi + m 2 ) 1 / 3 P 2 / b 3 

2 -0.462( m2 )V3 



m 1 + m 2 



(3) 
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where M_i is the mass loss rate of the secondary in units of 1CT 10 M yr _1 , n 2 o is the 
electron column density in units of 10 20 cm -2 , v$ is the speed in units of 10 8 cm s -1 , and 
-Porbh is the orbital period in hours. We adopt the parameters given above, the derived 
uh = 6.4 x 10 18 cm ~ 2 and 7 = 1. We finally choose v$ = 1, typical of winds driven by 
the radiation of a pulsar (Tavani & Brookshaw 1991) and ~3 times the escape velocity from 
the companion. We get M_io < 1.5. Even considering the uncertainty on this estimate, 
winds induced by the pulsar radiation are more typically 2-3 orders of magnitude weaker 
(Tavani & Brookshaw, 1991). On the other hand, if such a wind were independent from the 
association with the pulsar in a binary, it would be a common feature in the evolution of 
the single stars in NGC 6397. But the survival of lithium at the stellar surface is consistent 
only with rates smaller than 10~ 12 M Q yr" 1 (Vauclair and Charbonnel 1995) as a stronger 
wind, acting over the 10 Gyr of the cluster lifetime, would expose the stellar layers in which 
lithium has been destroyed by nuclear burning. The turnoff stars in NGC 6397 do show 
lithium at a level about normal for population II stars (Pasquini and Molaro 1996), putting 
strong upper limits to the mass loss rates of main sequence winds. 

Therefore, the mass loss rates which we derive are much more consistent with Roche 
lobe overflow driven by nuclear evolution of the secondary and orbital angular momentum 
mass loss, than with a possible wind from the secondary, suggested by D'Amico et al. (2001) 
and Ferraro et al. (2001). As our model requires that the mass lost is swept out by the 
pulsar radiation, the value of v$ given above is consistent with our hypothesis. 

3. Binary evolution in the Globular Cluster NGC 6397 

3.1. The HR diagram 

In order to choose coherently the input parameters for the possible evolution leading 
to PSR J 1740-5340 , it is important to take into account what we know of the general 
properties of the host cluster, which has been carefully studied down to the low end of the 
main sequence. Figure 1 shows the composite HR diagram of NGC 6397 in the plane M v 
versus V — I. We have plotted the sample by King et al. (1998) including the low mass 
main sequence (MS) and the white dwarfs (WDs), together with the sample by Cool (1997) 
for the turnoff and giant stars. The open circles identify the objects examined by Taylor et 
al. (2001) in the core of this cluster, to select objects which have been probably subject to 
binary evolution. In fact they recognize a helium WD sequence on the left of the MS, and 
many BY Dra stars on its right. The offset from the MS suggests indeed that these latter 
stars have suffered binary evolution. One of the Taylor et al. (2001) BY Dra candidates, 
plotted as a full dot in our figure 1, is indeed the optical counterpart of PSR J1740-5340 , as 
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discovered by Ferraro et al. (2001). On the observational HR diagram we show an isochrone 
of 12 Gyr for metallicity in mass fraction Z=0.006 and helium mass fraction Y=0.23 from 
Silvestri et al. (1998), complemented by the models by Baraffe et al. (1997) for masses 
< 0.5 M , while the track on the cluster white dwarfs is the 0.5 M Carbon Oxygen WD 
evolution by Wood (1995) as used in Richer et al. (1997). The isochrone of 12 Gyr which fits 
the cluster HR diagram implies that a mass of ~ 0.81 M is evolving at the cluster turnoff 
(TO), and that its TO luminosity is ~ 2.24 L . Different interpretations of the HR diagram 
morphology, assuming that the distance of the cluster is smaller, may lead to values of the 
TO luminosity down ~ 1.8 L Q . 



3.2. Choice of the initial parameters 

In order to choose the initial system parameters which leads to the present stage of 
PSR J1740-5340 , we make the following considerations: 

1) the location of the optical counterpart is at a luminosity similar to the TO luminosity (1.8- 
2.3 L , based on the previous discussion), but cooler than the TO, that is, at a radius larger 
than the TO radii. The pulsar luminosity is in fact estimated to be L PSR = 1.29 x 10 35 erg 
s _1 (see below), of which, in the hypothesis of isotropic emission, a fraction i? 2 ./ (2a) 2 impacts 
at the mass losing star surface (i? 2 and a being the radius of the secondary and the orbital 
separation respectively). Even in the most favorable case of a secondary filling its Roche 
lobe the fraction of the pulsar luminosity impacting the secondary is / = i? 2 lL2 /(2a) 2 = 
5.34 x 10 _2 [m 2 /(m 1 + m 2 )] 2 / 3 , where -Rrl2 is the Roche lobe radius of the secondary. We 
considered secondary masses ranging from 0.25 up to 0.81 M Q and corresponding NS masses 
from 2.0 down to 1.4 M (conservative mass transfer). We obtained values of /xi PSR ranging 
from 0.4 to 0.9 Irrespectively. Therefore the fraction of the pulsar luminosity reprocessed 
by the secondary is not very important for determining the HR diagram location of the star. 
The orbital period is moderately large, so any stage of mass transfer must have begun while 
the mass losing component was still not much evolved. 

2) We exclude the possibility that the star is a normal MS star losing a stellar wind at rates 
~ 1O _1O M /yr, for the reasons given at the end of Sect. 2. 

3) The pulsar spin up to 3.65 ms must be due to mass exchange in a previous evolutionary 
phase. Indeed the possibility that the NS was born with such a period seems unlikely as the 
spin-down age of PSR J1740-5340 is r = P spin /2P spin ~ 350 Myr (D' Amico et al. 2001), 
which is much shorter than the age of the globular cluster NGC 6397 in which PSR J1740- 
5340 is located. The easiest way is to attribute it to the previous evolution of the system. 
Therefore we model the initial parameter of the system as starting with a 0.85 M component 
and a NS component of 1.4M . NGC 6397 is a post core collapse globular cluster (Diorgovski 



- 8- 



& King 1986), in which a significant number of interacting neutron star plus main sequence 
binaries may have formed by tidal capture or exchange collisions (e.g. Davies, Benz and Hill 
1992, Di Stefano and Rappaport 1993). Implicitly, we are making the assumption that the 
0.85M Q component has been captured by the NS at such a separation to be able to begin 
mass transfer when it is not yet evolved as a giant. This is in line with what has been recently 
suggested by Podsiadlowski, Rappaport, and Pfahl (2002 CHECK!) that suggest, supported 
by the numerous ultracompact binaries found in globular clusters, that tidal capture seem 
to be a more likely way to produce this kind of system. 



3.3. The binary evolution 

We follow the binary evolution with the ATON1.2 code (D'Antona, Mazzitelli, & Ritter 
1989). The mass loss rate is computed following the formulation by Ritter (1988), as an 
exponential function of the distance of the stellar radius to the Roche lobe, in units of the 
pressure scale height. This method also allows to compute the first phases of mass transfer, 
during which the rate reaches values which can be much larger than the stationary values, 
due to the thermal response of the star to mass loss. The evolution of the system also 
includes orbital angular momentum losses through magnetic braking, in the Verbunt and 
Zwaan (1981) formulation, in which the braking parameter is set to / = 1. We also tested 
a case in which / = 2. 

We study four cases of evolution (see Table 1). In all cases, we assume that the secondary 
initial mass is M 2 = O.85M , the primary neutron star has Mi = 1.4M , and the orbital 
initial period is 14.27hr. The nuclear age of the O.85M is ~ lOGyr when it begins the 
Roche lobe overflow, so that it is only slightly evolved, and hydrogen is almost exhausted in 
its (radiative) core. In case 1, we assume a 'standard' conservative evolution. This produces 
mass transfer rates of ~ 2 — 4 x 10~ 10 M yr -1 . The system would thus appear as a low mass 
X-ray binary (LMXB). The X-ray emission expected from the NS is ~ 2 x 10 36 erg/s. The 
evolution is typical for such systems (e.g. Bhattacharya & van den Heuvel, 1991). The end 
product will be a very low mass helium white dwarf (M 2 = O.246M ) in a much wider orbit 
(-Porb = 119 hr) with a 2M neutron star. Being the hydrogen not yet completely exhausted 
in the stellar core when mass transfer begins, the initial binary period we have chosen is 
indeed very close to the "bifurcation" period below which the orbital evolution proceeds 
towards shorter binary periods. Note that in our case the bifurcation period is shorter than 
that quoted by Podsiadlowski, Rappaport, and Pfahl (2002) for a binary composed by a 1.4 
M NS plus a 1.0 M secondary, namely P or bbif ~ 18 hr. This is due to the fact that the 
main sequence radius of our 0.85 M population II secondary is ~ 20% smaller than the 
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radius of a population I, 1.0 M companion for comparable consumption of the hydrogen 
in the core. Since for a Roche-lobe filling secondary we have that AP or b/P rb — 2.2AR2/R2, 
the bifurcation period for the population I, 1.0 M companion will be ~ 44% longer than 
in our case and therefore consistent with the result of Podsiadlowski, Rappaport, and Pfahl 
(2002). 

Having in mind the oberved system PSR J1740-5340 , we followed another case of 
evolution (case 2) in which we assumed conservative evolution up to an orbital period of 
~ 26 hr, and then we assumed that the matter lost from the secondary is all lost from the 
system at the inner lagrangian point LI with its specific angular momentum. This alters 
slightly the conservative evolution, leading to marginally different mass loss rates and final 
parameters for the system when all the seconday hydrogen envelope is lost and the star 
becomes a white dwarf (see Table 1). 

When mass accretion on the NS is assumed, and the binary evolution suffers a LMXB 
phase, it is also important to consider how the X-ray phase would affect the binary evolution 
we are considering, as a fraction R\/(2a) 2 < 2.7 x 10~ 2 of the X-ray luminosity (see point 
1 above) impacts at the secondary surface. There have been many attempts to model the 
effect of this irradiation on the binary evolution, when it has been realized that spherically 
symmetric illumination in LMXBs can have a dramatic effect on the structure of low mass 
secondaries having convective envelopes. In fact, X-rays block the intrinsic stellar flux 
produced by the nuclear reactions, leading to a stellar expansion, which affects the mass 
transfer rates (Podsiadlowski, 1991; Harpaz & Rappaport, 1991). In reality, illumination 
will be non-spherical, because only one hemisphere is affected - in fact the donor star is 
tidally locked as it is filling its Roche-lobe during this evolutionary phase -, and its effect 
will depend on the depth at which the energy is deposited below the photosphere and how 
fast it can be transported to the cool side of the star. If the circulation time is shorter 
than the cooling time of the heated matter, we go back to the spherical illumination case. 
Otherwhise, the global secular evolution is not much affected, but the mass transfer occurs 
in "outbursts" whose duration and peak mass transfer rates depend on the efficiency of 
circulation (e.g. Hameury et al. 1993, Harpaz and Rappaport 1994, Vilhu et al. 1994, 
D'Antona 1994). All the existing computations refer to the evolution of LMXBs with main 
sequence companions. We are now dealing with a case of evolution in which the mass losing 
component is more evolved, and in addition has no important convective envelope at the 
time the mass transfer begins, so that it is important to have an idea of the possible effects 
of X-ray illumination in our specific case. 

On the other hand, self-consistent modelling is really very difficult, as the feedbacks are 
not easy to describe both physically and numerically. However we follow the evolution of an 
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ideal limiting case, considering a system in which the secondary is affected by a fixed value 
of heating luminosity Lh = 5 x 10 34 erg s -1 (case 3). We include the effect of illumination 
as described in D'Antona and Ergma (1993), following Tout et al. (1989). In this simplified 
schematization, the star is immersed in the X-ray radiation bath, and the total luminosity 
L to t which it must radiate is the sum of the stellar luminosity plus the heating luminosity 
Lh, from which an 'irradiation temperature' T irr = (Lh/^naRl) 1 ^ 2 is defined. The stellar 
luminosity and T e g then are related by: 

L = L tot -U = AnaRKT^ - 7*.) (4) 

Consequently, the stellar T cS becomes hotter due to the irradiation, and the star rapidly 
evolves in the HR diagram at a location determined by the amount of irradiation allowed. 
The phases of mass transfer are slightly altered by the new system conditions, but are globally 
very similar to case 1. However, as the star loses mass, its radius becomes larger than the 
radius of the standard sequence. This difference amounts to ~ 20% at the orbital period 
of PSR J1740-5340 . Of course the X-ray irradiation is not self consistently described by 
this model. In fact, any fluctuation in the mass transfer rate will be amplified. If the mass 
transfer decreases further, irradiation decreases too, the radius also will decrease, trying to 
reach its non irradiated equilibrium value and the mass transfer rate decreases more. We 
can foresee that the LMXB phase in fact can alternate phases of mass transfer and detached 
phases, as it was predicted for the LMXB having MS secondaries (Hameury et al. 1993, 
D'Antona 1994, Harpaz and Rappaport 1995). 

To test the sensitivity of the mass loss rate to the assumptions made concerning magnetic 
braking, we follow another evolution (case 4) in which / = 2 and the mass transfer is non 
conservative for P orb > 17.36 hr. 

3.4. The final phases of evolution 

Figure 1 shows the evolutionary paths of case 2 and 3 in the HR diagram. The theoretical 
values of luminosity and T c g- are converted into the observed magnitudes V and I by means 
described in Bessell, Castelli and Plez (1998). The irradiated track (case 3) represents the 
evolution of the accreting progenitor of PSR Jl 740-5340 and therefore does not reproduce 
the location of PSR J1740-5340 for which the accretion is inhibited by the pulsar radiation. 
The track passing through the optical counterpart of PSR Jl 740-5340 is the case 2 evolution. 
The track corresponding to case 1 (purely conservative evoluton) is very similar, and is not 
shown. 
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All the sequences are evolved until the mass loss phase finally ends with the stellar 
remnant evolving into the white dwarf region as low mass helium white dwarfs. We see that 
the most luminous three objects among those identified by Taylor et al. (2001) as helium 
WDs actually may be the end-products of such an evolution. Notice however that further 
study is necessary to assess whether the lower luminosity objects are helium white dwarfs. 
In fact, unless their colors are peculiar, they seem too cool to have the radius expected, even 
for the lowest possible helium WD masses remnant of mass exchange evolution, which are 
actually ~ 0.2 M (D'Antona et al. 2002, in preparation). The open triangles show the 
location at which sequence 3 achieves a total age of 17, 20 and 26 Gyr (~ 10.5 Gyr of which 
have been spent in the phases previous to mass exchange). The point at 20 Gyr, which 
corresponds more closely to the luminosity of the lowest Taylor et al. objects, is at T e g= 
8900 K. The full squares along sequence 2 correspond to ages from 13 to 20Gyr in steps of 
lGyr, and a last square indicates the location at 25Gyr. We also notice that the cooling 
of such objects (also in the present models) is dominated by proton-proton burning of the 
remnant hydrogen layer (Driebe et al. 1998, Sarna et al. 1998, Schonberner et al. 2000), 
and that the low initial metallicity of the system prevents hydrogen shell flashes. 



4. The evolution of the progenitor system of PSR J1740-5340 

The computation of case 3 has shown that the binary evolution which we are describing 
is not dramatically altered by the irradiation due to the LMXB phase: in particular, the final 
evolution is very similar, although the white dwarf remnant mass is slightly smaller. However, 
it helps to predict that the LMXB phases can be alternated with phases in which the system 
remains detached, as already had been suggested for systems having MS companions. Let us 
now consider the evolution of the system, taking also into account the pulsar's behavior. Any 
time the system detaches, the radio pulsar will switch on, but it will be quenched again when 
the mass transfer resumes. In the HR diagram the position of the secondary component will 
shift from its 'irradiated' position during the LMXB mass transfer phase (track 3 on the left 
of the MS) to its 'standard' position along track 1. However, as discussed in Burderi et al. 
(2001a,b), when the pulsar is spinning sufficiently fast and the orbital period is sufficiently 
long (i.e. the orbital separation is large), the radiation pressure exerted by the pulsar at the 
inner inner Lagrangian point is larger than the pressure exerted by the matter overflowing 
the Roche lobe even if the mass transfer rate recovers its secular value dictated by the nuclear 
evolution of the companion. 
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This condition is verified for orbital periods longer than 1 

/Wit = 0.75 a~ 54 / 25 n^f m 107/50 K + m^" 1 / 2 

g{mi,m 2 ) fi 2 e ^ sp m-3 M -io n W 

where a is the Shakura-Sunyaev viscosity parameter, n .6i5 = n/0.615 ~ 1 for a gas with 

solar abundances {n is the mean particle mass in units of the proton mass m p ), g (mi,m 2 ) = 
/ \ 1/3 

1-0.462 ( mi + 2 m J , /u 26 is the magnetic moment of the NS in units of 10 26 G cm 3 (fx = B S R 3 
with R and B s NS radius and surface magnetic field along the magnetic axis, respectively), 
and P sp in _3 is the NS spin period in milliseconds. If P or b > P or h crit the system will remain 
in the radio-ejection phase during all the subsequent binary evolution. In this case, the 
matter overflowing the Roche lobe will be accelerated by its interaction with the pulsar 
radiation and ejected by the system. If we can assume (but at least ballistic simulations 
should be done) that this matter leaves the system with the specific angular momentum at 
the inner Lagrangian point, I = (a — -Rrx2 — ?r) 2 x (27r/P or b) (r\ is the distance between 
the NS and the center of mass of the system) which in our case is quite close to the specific 
angular momentum of the system, the orbital evolution is indeed very similar to the more 
conventional evolution with conservative mass transfer to the primary (as confirmed by our 
numerical computations of the conservative case 1 and the non conservative case 2). 

The optical component of PSR J1740-5340 indeed seems compatible with the evolution 
we have suggested. It lies along case 2 evolution, that not including irradiation, as, in fact, 
the pulsar luminosity is not important as irradiation source for the binary. 



1 Equation 5 has been derived in Burderi ct al. (2001a, b) equating the pressure of a Shakura-Sunyaev 
accretion disc with the radiation pressure exerted by the pulsar. Once the disc has been swept away the 
pressure of the flow overflowing the inner Lagrangian point should be used instead of the pressure of the 
disc. However it is possible to demonstrate that these pressures are comparable, which supports the validity 
of equation 5. Let us suppose that the flow of matter overflowing the inner Lagrangian point has a cross 
sectional area Ar ow = 7r(<5_RriLi) 2 with 5 < 0.1. Moreover the inward speed of the flow is va ow — f3vs 
where vs is the free fall speed and j3 < 1 since the Coriolis force deviates the flow motion from the radial 
direction. The ram pressure of the flow is Pfi ow = pf[owv\ owl where pa ow is the density of the flow that is 
related by the continuity equation to the mass loss rate. The pressure of a Shakura-Sunyaev accretion disc 
is Pss — PsskT/ (nm p ), where pss is the density of the disc that is, also in this case, related by the continuity 
equation to the mass loss rate, k is the Boltzmann constant, and T is the temperature of the disc at its 
midplane (see e.g. Frank, King, and Rainc 1992). After some algebraic manipulation, using the standard 
relations for a Sakura-Sunyaev accretion disc, we get Pfl ow = -Pss x 100(/i/r)(4/v / 2)a(/3/(5g j), where h is the 
half thickness of the disc and So.i = 5/0.1. As a good approximation we can take h/r <~ 1.7 x 10~ 2 as the 
dependencies on the disc parameters are rather weak (see again Frank, King, and Raine 1992). Therefore 
we have P flow ~ 5a((3/5$ A )Pss- 
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5. Summary and conclusions 

We have considered the evolution of possible progenitors of the binary MSP PSR J1740- 
5340 in the Globular Cluster NGC 6397. We can reproduce the HR diagram location of the 
optical companion, starting mass transfer to the NS from a hypothetical secondary of mass 
0.85 M , slightly evolved off the MS when mass transfer begins. From equation 5, adopting 
a — 1, no. 615 = 1, rn\ = 1-8, m 2 = 0.45, U26 = 7.7, P sp m -3 = 3.65, and M_i = 1, the critical 
period P or bcrit to reach the 'radio-ejection' phase is ~ 39 hr, not very different from ~ P or b, 
given the large uncertainties in the M_io estimate. Thus, PSR J1740-5340 is possibly in the 
radio-ejection phase. Given that P or bcrit oc P S pf n u~ 4 - 8 M 2 - 04 the fact that P or bcrit ~ Porb is 
compelling. 

In conclusion we suggest that: 

i) orbital evolution calculations shows that a slightly evolved 0.85 M secondary orbiting 
a NS can transfer mass to the NS, and reaches a stage in which its mass is reduced to 
~ 0.45 M , and its optical location in the HR diagram is then compatible with the 
recently detected optical counterpart of PSR J1740-5340 ; 

ii) PSR J1740-5340 might represent a system whose evolution has been envisioned by 
Burderi et al. (2001): the spin and the magnetic moment of the pulsar may keep 
the system in a radio-ejection phase in which accretion is inhibited by the radiation 
pressure exerted by the pulsar on the overflowing matter while the mechanism that 
drives the Roche lobe overflow from the companion is still active, thus causing an 
intense wind which would be very difficult to explain otherwise. This evolution seems 
to be the only viable possibility to explain the long lasting eclipses and the strong 
intensity variation randomly occurring in the radio emission. 

As a final remark we note that P or b ~ Porbcrit suggest the interesting possibility that 
this system could swiftly switch from the present radio pulsar phase to an accretion phase 
in which it should be visible as a Lx ~ 10 36 ergs s -1 LMXB. 
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enlightening discussions during the preparation of this work. This work was partially sup- 
ported by a grant from the Italian Ministry of University and Research (Cofin-99-02-02). We 
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V-I 

Fig. 1. — HR diagram of the Globular Cluster NGC 6397, with an isochrone superimposed 
(see the text for details). Starting close to the turnoff at the open square, we show a standard 
evolution (case 2) of the companion of the NS, having an initial mass of O.85M and an initial 
orbital period of 14.4 hrs. (track evolving first towards the right of the figure, and passing 
through the observational counterpart of PSR J1740-5340 - full dot), and the evolution of 
case 3, which includes a fixed irradiation luminosity (track to the left). Both sequences 
end into the WD evolution of a 0.246 M Q (case 2) and 0.216 M (case 3). Both sequences 
are followed along the WD cooling, which is dominated by residual proton proton burning 
lasting for more than a Hubble time. Ages of 17, 20 and 26 Gyr are labelled as triangles 
along the cooling track of case 3. Ages from 13 to 20 vin steps of 1 Gyr, plus a last point at 
25 Gyr, are labelled as squares for case 2 
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Table 1. Computed evolutions: M lin =1.4M , M 2in =O.85M , P in =14.72hr 



Case 


M NS fin 


M WD 


Porb 


AML 


Mass transfer 


u 


log M av 




M 


M 


hr 




modalities 


erg/s 


M /yr 


1 


2.004 


.246 


119 


MB f=l 


Ml+M2=c 





-9.5 


2 


1.72 


.247 


120 


MB f=l 


Ml=c from 





-9.5 












P=25.98hr 






3 


2.03 


.217 


278 


MB f=l 


.1/1 • .1/2 C 


5xl0 34 


-9.9 


4 


1.52 


.261 


193 


MB f=2 


Ml=c from 





-9.3 



P=17.36hr 



